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Abstract — Shape is one of the most important visual attributes to characterize ob- 
jects, playing a important role in pattern recognition. There are various approaches to 
extract relevant information of a shape. An approach widely used in shape analysis 
is the complexity, and Fractal Dimension and Multi-Scale Fractal Dimension are both 
well-known methodologies to estimate it. This papers presents a comparative study 
between Fractal Dimension and Multi-Scale Fractal Dimension in a shape analysis 
context. Through experimental comparison using a shape database previously classi- 
fied, both methods are compared. Different parameters configuration of each method 
are considered and a discussion about the results of each method is also presented. 
Keywords — Complexity; Shape Analysis; Fractal Dimension; Multi-Scale Fractal 
Dimension; Fourier Transform. 



1 Introduction 

In pattern recognition and image analysis, shape is one of the most important visual 
attributes to characterize objects. It provides the most relevant information about an 
object in order to perform its identification and classification. It is a classical problem, 
and literature presents a large amount of techniques to extract information related to 
shap e geometric aspect, allowing to separate and to label different parts of an image 
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An approach widely used in shape analysis applications is the study of shape through 
its complexity. In this analysis, complexity is straight related to the irregularity pattern 
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presented by the shape under analysis and, respectively, the amount of the space the 
shape occupies m. 

An interesting way to estimate the complexity of an object is using the Fractal 
Dimension il2ll . Different of topological dimension, which is an integer number, the 
Fractal Dimension is a fractionary value that describes how irregular is an object and 
how much of the space it occupies. Bouligand-Minkowski method is one of the most 
accurate methods to compute Fractal Dimension. Is is based on the study of the shape 
influence area computed by shape dilation 

An alternative method for estimating shape complexity is the Multi-Scale Fractal 
Dimension. It consists in estimating a curve that represents the changes in complexity 
shape as we change the visualization scale. Different of Fractal Dimension, which is a 
numeric value, this appro ach produces a curve which performs a more accurate shape 
discrimination |3 J 3 131. 



In this paper, we propose to evaluate both Fractal Dimension and Multi-Scale Frac- 
tal Dimension in a shape analysis context. For this, a shape database is used and ex- 
periments are performed considering different configurations of each method. This 
paper is organized as follows: Section|2]describes the Fractal Dimension method while 
section [3] shows how to compute the Multi-Scale Fractal Dimension. In section]?] is 
showed how to use the Multi-Scale Fractal Dimension as a shape signature. The sig- 
nature making process is based on Fourier analysis and it is described in section ]5] 
Section]6]describes the methodology adopted by this work. Results and conclusion are 
finally given in sections]?] and ]8] respectively. 



2 Fractal Dimension 

Fractal Dimension is a measure of how fragmented a fractal object is, and it may be 
understood as a characterization of its self-similarity It is a non-integer number 
that quantifies the density of fractals in the metric space and it is a way to identify how 
complex a fractal is, in order to compare it with another. Fractal geometry has various 
approaches to compute the Fractal Dimension of an object. These approaches can be 
classified as belonging to the Hausdorff-Besicovitch Dimension (like the BoxCounting 
and Dividers methods) or to the Bouligand-Minkowski Dimension (Minkowski Frac- 
tal Dimension method), where the last is the one which produces the most accurate 
and consistent results for Fractal Dimension 1 3. 3 0] . Bouligand-Minkowski Fractal 



Dimension method is based on the study of the influence area created by shape di- 
lation using a disc of radius r (Figure JT]). Small modifications on the shape produce 
modifications in the computed influence area |T3, Hf]. 



Consider A E the shape under analysis, the dilation of A, A{r), is defined as 
the set of points in such distance from A is smaller than or equal to r: 



A{r) ^ {x e R^\3y <E A : \x ~ y\ < r} ^ 
This dilation can also be defined as: 

A{r) = U Br{x), 

xeA 
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Figure 1 : Dilation of a shape with a disc of radius r. 



where Br{x) is a disc of radius r. The influence area, A{r), and the radius, r, 
follows the relation: 

A{r) =^r2-^. 
So, the Fractal Dimension can be estimated as: 



D = 2 - lim 



log A{r) 



>o log r 

where D is the Fractal Dimension estimated by the Bouligand-Minkowski method. 
Through line regression of log-log curve A{r), it is possible to calculate a line with a 
slope, where D = 2 — a is the Fractal Dimension of the shape using the Bouligand- 
Minkowski method llT6ll . 

Computing the influence area of an object is a task of high processing cost. One 
possibility to optimize this task is using the Euclidean Distance Transform (EDT), 
which attributes to the pixels of a binary image the minimum distance among pixels 
from the image object to the pixels from the image background. The distance func- 
tion considered duringthe EDT computing is the Euclidean distance, due its rotation 
invariance (Figure |2|l 1411411. 




Figure 2: Euclidean Distance Transform (EDT) applied over an binary image. 
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3 Multi-Scale Fractal Dimension 



The main problem with Fractal Dimension in nature shape characterization is that they 
are not fractals (or self-similar). All objects have a finite size, and it implies that their 
complexities go to zero as long as observation scale increases. An interesting inter- 
pretation of the fractal behavior of a shape is how the dilation, using the Bouligand- 
Minkowski method, occurs in different points of the shape. Some shapes allow their 
poin ts to be freely dilated while in other points this dilation is saturated in some radius 

Mi. 

This behavior gives to the Bouligand-Minkowski log-log curve a richness of details 
that can not be expressed by just a numeric value, as performed using line regression. 
Using the derivate, it is possible to find a function that binds the Fractal Dimension 
changes to the dilation radius changes (Figure O 10, [13, [Hi ■ This function is called 
Multi-Scale Fractal Dimension (Af FD), and it is defined as: 

du(t) 

MFD = 2 ^ 

dt 

where du{t)/dt is the derivative of log-log curve u{t) calculated by Bouligand- 
Minkowski method. 

11 I . . . . 1 1.4 I . . . . . — 
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Figure 3: (a) log-log curve from Bouligand-Minkowski method, (b) Multi-scale Fractal 
Dimension. 



In order to compute the MFD, it is necessary to calculate the derivative of u{t). 
This is performed using the derivative property of the Fourier Transform. This property 
allows to compute the derivative curve in the spectrum and it has a better performance 
when compared with numeric methods, once it considers all data points during the 
derivative computing 10, [3 • An important detail that requires attention is that deriva- 
tive methods have a tendency to emphasize high frequency noise. So, it is necessary to 
use a low pass filter, like Gaussian filter, in order to reduced the this noise 10. [I]- The 
derivative based on Fourier Transform can be expressed as: 

^ = F-^ {F {u{t)} F {g,{t)} {j2nf )} 
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with 



1 f-t^ 



where t and u{t) are, respectively, the logarithm of radius and the influences area 
from Bouligand-Minkowski method, / is the frequency, j is the imaginary number and 
ga{i) is the Gaussian function with standard deviation a [9, .l(\ll3ll . 

Although the advantages of the derivative property of the Fourier Transform, some 
important aspects must be considered when computing the MFD curve. These as- 
pects regards to the log-log curve u{t) behavior and they are necessary to achieve a 
good shape discrimination. At first, it is necessary to provide a curve with a good 
sampling and uniform interval. This is performed when the initial points of the curve 
are not considered, once they present low sampling (Figure |4|i, followed by a linear 
interpolation, which is performed by filling the space between each two points of the 
sampled curve by its average point. 
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Figure 4: In light gray: low sampling region where the points are not considered. 

Another problem present in the derivative property of the Fourier Transform is the 
discontinuity of the method in the curve limits (Figure|5]l. This phenomenon is known 
in literature as Gibbs phenomenon jst]. This phenomenon is due to the fact of the 
Fourier transform do not converge uniformly in discontinuities. 

An effective solution for this problem is to use a scheme of duplication and reflex- 
ion of the curve, so that, it is possible to make it continuous (Figure |6]l. This scheme 
provides a continuous curve for the interval [2iV — 1, 3iV], where N is the length of the 
original curve. 
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4 Multi-Scale Fractal Dimension as a shape signature 



Image signature is defined as a simplified function or matrix that is able to represent 
or characterize the original image. In general, techniques used for this task apply a 
T : /•^ — > i? or T : — > / transformation, where the computed feature vector 
represents the original image in a simplified way. More theory and examples can be 
found in flSll. 

This paper presents a study about the possibility of using the MFD technique 
as a shape complexity signature. Different of true fractals, images from nature are 
not self-similar They may look different whether we change the visualization scale 
and, consequently, their Fractal Dimension is also dependent on the used scale. The 
MFD technique allows to study the shape behavior through scales. As a result, a 1-D 
function is calculated for a given shape, where this function describes how the Fractal 
Dimension of that shape changes as we change the scale (Figure]?]). 

1.6 I 1 1 1 1 1 1 1 
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Figure 7: MFD curve as a shape signature. 



5 Fourier Descriptors 

The Fourier transform is a useful technique that allows to study the behavior of a signal 
through its frequency spectrum. It is widely used in pattern recognition tasks, once it 
presents a great number of advantages, such: noise tolerance, analysis of a signal into 
different groups of frequencies and easy data normalization, yielding data which is 
invariant to rotation, translation and scale (HQl- 

When it is applied over a signal u{t), the Fourier transform obtains the complex 
components U{f): 
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U{f) = / u{t)e-'^-^dt, 

which are the signal representation in the frequency spectrum. This approach al- 
lows to split the original signal into different groups of frequencies, each one with 
distinct signal features. 

Lower frequency coefficients, for example, are associated to the portion of the spec- 
trum that describes the most relevant information about the signal behavior, whereas 
higher frequency coefficients possess information about noise and details presented in 
the signal [TO. For the MFD, only the lower frequency coefficients are considered. 
The Fourier descriptors, DU{f), are computed from the magnitude of the selected low 
frequency coefficients: 

DU{f) = \\U{f)\\. 

A normalization is also performed in such descriptors, and it is done as follows: 

DU{f)=DU{f)/DU{l). 

This normalization process increases the tolerance of the descriptors to distur- 
bances in the original signal, such as changes in scale, translation and rotation. 

6 Experiment 

Considering the possibility of using the MFD as a shape signature, an experiment 
based on image classification is performed to evaluate its efficiency. 

An image database with 1352 images previously classified into 26 classes of 52 
images is employed. Each class corresponds to an uppercase letter of the occidental 
alphabet (26 letters latin alphabet). Images from each class are grouped into 4 groups 
of 13 images, where each group corresponds to a different random noise level. This 
allows to evaluate the noise tolerance of the method. 

In order to achieve the best performance of the method, various parameters values 
have been tested. These parameters refers to the dilation radius r used in Bouligand- 
Minkowski method and the smoothing level a used in the MFD calculation. For each 
different configuration, a new MFD curve is computed for each sample in the image 
database. An average MFD curve is also computed for each class using 12 samples 
of it (3 samples for each noise level), and these used curves are discarded in order to 
not affect the classification results. So, the test database is compound of 1040 image, 
divided into 26 classes of 40 images. Each class is also divided into 4 noise levels, each 
one with 10 images. Classification process is performed by Euclidean distance, where 
it is verified the distance among the remainders curves and the average M FD. 

7 Results and Discussion 

In this section some properties of MFD curve are discussed, as well as its performance 
according to its parameter. Besides, we also compare the results from MFD and 
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Figure 8: Example of noise levels: (a) Original Shape; (b) Level 1; (c) Level 2; (d) 
Level 3; (e) Level 4. 



Fractal Dimension, and the advantages and disadvantages of using Fourier Descriptors 
from MFD curve for shape classification. 

7.1 Method parameters analysis 

In order to attain the best method performance, various parameters values have been 
tested. The tested parameters values are cr e {10,15,20,25} and 
r e {10,15,20,25,30,35,40,45,50,75,100,125,150,175,200,225}. The results 
for each different configuration showed the characteristics and properties present in 
MFD curve. 

We note an increase in success rate as r increases (Figure |9]l. However, the rela- 
tion between success rate and the dilation radius presents a log-curve behavior So, 
from an specific radius value, the success rate becomes constant. The explanation for 
this behavior lies in the fact that as r increases, the dilation produced by Bouligand- 
Minkowski method makes the shape aspect becomes more similar to a point, whose 
Fractal Dimension is equal to zero. This behavior is represent in MFD curve as its 
tendency to zero according to r increase. From an specific r, none relevant information 
about the shape is added to MFD curve, except for its zero-tendency, what keeps the 
success rate stable. We also note that this behavior do not depends on <j value (Figure 



9 



Besides r parameter, which refers to dilation radius, we realize that cr value also 
carry out some influence over the success rate. The parameter a refers to the smoothing 
level of MFD curve and it is responsible for the presence of more or less details in 
the curve (Figure [TTTi. We note higher values for this parameter produce a excessive 
smoothing and it suppresses important details of the curve. As a consequence, the 
success rate decreases, independent of the value of parameter r (Figure [T2]). 

In this experiment, the best performance is achieved when used r = 100 and a = 

10. 

100 I 1 1 1 1 1 
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Figure 9: MFD success rate as radius increases. 



7.2 Fractal Dimension versus Multi-Scale Fractal Dimension 

As previously discussed, the main problem with Fractal Dimension in shape character- 
ization is that it does not work properly with real world restrictions. As we increase 
the visualization scale, the Fractal Dimension of the object goes to zero. In the other 
hand, MFD curve emphasizes the details present in log-log curve computed from 
Bouligand-Minkowski method. It allows to produce a more detailed representation of 
how complexity changes according to the scale. This representation is much more con- 
sistent for analysis and less dependent on noise interferences than a simple numeric 
value as calculated by u(t) line regression. 

Figure [13] shows the classification results for both methods. Fractal Dimension 
(FD) and Multi-Scale Fractal Dimension (MFD). For this experiment, various ra- 
dius values have been tested and, in MFD case, an = 10 was considered. 

We note a higher success rate for MFD when compared with FD, independent 
of radius values used. In FD method, only an numeric value is calculated from u{t) 
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Figure 10: Example of MFD computed for different dilation radius r, using a = 10: 
(a) r = 75; (b) r = 100; (c) r = 125; (d) r = 150. 

and used for shape characterization. This becomes the method more sensitive to small 
variations or noise in the shape. The results from DF are also more unstable than 
MFD. The explanation for this lies in the zero-tendency presented in the Bouligand- 
Minkowski method. As the dilation radius increases the Fractal Dimension of any 
object goes to zero, what becomes the classification process more difficult and inaccu- 
rate. 

Table s[T] and |2] show, respectively, the confusions matrix for DFM and DF when 
considered all images in the test database and r = 100. 

7.3 Noise Tolerance 

An important characteristic to be evaluated in any method is its capacity to work prop- 
erly with samples that present sort of distortion or noise. In general, images under anal- 
ysis may present distortions from measuring errors (measures from physical world), or 
noise added to original data (e.g., interference added during data transmission), that 
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Figure 12: MFD success rate for different a values. 
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Table 1: Confusion Matrix for Multi-Scale Fractal Dimension method, considering 
r = 100 and a = 10. Success rate = 94.62 %. 
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11 3 
1 7 



18 2 
1 2 



8 10 
4 4 



Table 2: Confusion Matrix for Fractal Dimension method, considering r = 100. Suc- 
cess rate = 47.40 %. 



difficults the analysis and classification processes. 

By analysis of the results from M FD curve, we note this method present a good 
noise tolerance. This tolerance comes from using a Gaussian low filter during the 
computing of MFD curve. The Gaussian filter acts reducing the importance of the 
informations in the high frequency region of the spectrum. Once noise is a typical high 
frequency information, it is removed from the original data through this data filtering 
process. We also note this filtering process presents a better performance on images 
that present a intermediate noise level, i.e., levels 2 and 3 of the used noise (Figures 
[8]: and [8]l, respectively) . In other hand, we also realize a stressed decrease in the 
performance of the method as the noise exceed these noise levels. High noise levels 
change the geometric patterns in the shape, what difficults the classification process 
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Figure 13: MFD and FD success rates as dilation radius increases, 
and, as a consequence, decreases the success rate of the method (Figure [l4]l. 
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Figure 14: MFD success rate at different noise levels and dilation radius. 



7.4 Analysis using Fourier descriptors 

As previously discussed, the main advantage of MFD curve in pattern recognition 
is its capacity to represent an object by a curve that binds its complexity along the 
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visualization scale. Nevertheless, the M FD curve presents a high cost to be evaluated 
once its number of descriptors increase as the radius r increases (Figure [TS]). Besides 
that, not all information in the curve is relevant for shape classification. 
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Figure 15: Number of descriptors in the M FD curve for the radius used. 

In order to simplify the analysis process the Fourier transform is applied over 
MFD curve, so descriptors that represent the main behavior of the curve are com- 
puted, as mentioned in Section|5] 

One important point to be considered is the number of descriptors to be computed. 
The number of descriptors have been defined through analysis of the Euclidean distance 
between the original M FD curve and its respective set of descriptors. Figure[T6l shows 
the variation of this distance as we increase the number of Fourier descriptors. We 
realize as we increase the number of descriptors smaller becomes the distance between 
the MFD curve and its descriptors. However, this distance becomes stable when the 
number of descriptors is higher than 50. This shows that most of the main information 
lies in the 50 first descriptors. 

Once it was defined the number of descriptors that best represent the curve behav- 
ior, the performance of these descriptors was compared with the original curve in the 
shape classification experiment. For this, we considered the MFD curve calculated 
for different dilation radius, r, and a — 10. Figure [17] shows the success rate for both 
Fourier descriptors and MFD curve. 

Using Fourier descriptors is a way to simplify the analysis step, once it reduces 
the amount of information to be analyzed. Otherwise, suppressing part of curve in- 
formation, specially high frequency information, produces a smoothing effect over the 
MFD curve, which decreases the distance between classes and, consequently, de- 
crease the success rate. In spite of the large amount of information discarded during 
the process of Fourier descriptors computing, the decrease in success rate is minimum 
for r > 50. This range of r values is where the method has presented its best results. 
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Figure 16: Distance between MFD and its Fourier descriptors as the number of com- 
puted descriptors increase. 
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Figure 17: MFD and its Fourier descriptors success rate. 

and it becomes evident the great capacity of MFD curve, in association with Fourier 
descriptors, to shape classification. 



8 Conclusion 

This paper have presented an experimental comparison between Fractal Dimen- 
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sion and Multi-Scale Fractal Dimension in a shape analysis context. Fourier descrip- 
tors have also been performed over Multi-Scale Fractal Dimension curve in order to 
simplify the analysis step. A previously classified shape database was employed and 
different configurations of each method have been tested. 

Results have shown that Multi-Scale Fractal Dimension performs a more accurate 
discrimination of the shapes. It also presents a "good" noise tolerance, what is corrob- 
orated by experimental results. Experiments also have shown that Fourier descriptors 
computed from Multi-Scale curve hold the main information of the curve, with minimal 
loss of information, while it reduces drastically the number of descriptors necessary to 
shape classification. 
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